Fate of the Higgs mode near quantum criticality 



Snir Gazit, Daniel Podolsky, and Assa Auerbach 
Physics Department, Techmon, 32000 Haifa, Israel 
(Dated: December 18, 2012) 

We study the relativistic O(N) field theory near the quantum critical point in 2+1 dimensions 
for TV = 2 and TV = 3. The scalar susceptibility is evaluated by Monte Carlo simulation. We show 
that the spectrum contains a peak associated with the Higgs mode, which remains well-defined all 
the way to the critical point. The fidelity of this peak and the amplitude ratio between the critical 
energy scales on both sides of the transition are computed. 
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Spontaneously broken continuous symmetry in con- 
densed matter produces collective modes. In addition to 
Goldstone modes, an amplitude (Higgs) mode is some- 
times expected at finite energy [TJ [5] . Higgs oscillations 
have been measured in e.g., the superconductor NbSe2, 
[U 03 H] , the dimerized antiferromagnet TICUCI3 [5] , and 
charge density wave compounds [SHE]. 

In the absence of gauge fields, the massive Higgs mode 
decays into massless Goldstone modes, broadening its 
spectral line. For relativistic O(N) models in 3+1 dimen- 
sions, the Higgs mode becomes an increasingly sharper 
excitation the closer one gets to the quantum critical 
point (QCP) [5]. This is a consequence of the fact that 
the QCP itself is a Gaussian fixed point. In contrast, 
in 2+1 dimensions (d = 2) the QCP is strongly coupled 
|10j , and there is no a priori reason to expect the Higgs 
mode to survive near criticality. 

Recent interest in the fate of the Higgs mode in 2+1 
dimensions has led to new theoretical and experimental 
results. The visibility of the Higgs peak has been shown 
to be sensitive to the symmetry of the probe jTTJ IT2"] : 
The longitudinal susceptibility diverges at low frequen- 
cies as w _1 [T3"rtPo] . This is due to the direct excitation of 
Goldstone modes which can completely conceal the Higgs 
peak. In contrast, the scalar susceptibility [T2] rises as a; 3 
and its Higgs peak is much more visible, even at stronger 
coupling. 

Indeed, in recent experiments of cold bosons in an op- 
tical lattice [TB] , the Higgs mode has been detected in the 
scalar response, in the vicinity of the superfluid to Mott 
insulator transition. Further large N analysis [T7], and 
numerical simulations of the Bose Hubbard model [18] 
(N=2), have suggested that the Higgs peak is still visible 
as it softens toward the QCP. However the ultimate fate 
of this peak in the critical region demands simulations on 
much larger systems. 

The question to be answered is whether the Higgs 
mode is well defined arbitrarily close to criticality, in the 
sense that one can identify a peak in the scalar suscepti- 
bility even as its energy scale decreases to zero. 

In this Letter we answer this question in the affirma- 
tive. We compute the frequency dependent scalar suscep- 
tibility for the relativistic 0(2) and 0(3) models. Lorentz 
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FIG. 1: Critical energy scales near the quantum phase tran- 
sition in relativistic O(N) field theory for N — 2,3. Sg = 
(g — g c )/g c is the dimensionless tuning parameter, ran is the 
Higgs peak energy (mass) in the ordered phase Sg < 0, and A 
is the gap in the disordered phase Sg > 0. Solid lines describe 
the critical behavior m H = B-\5g\" N and A = B+\5g\" N . In- 
set: Phase diagram in the microscopic g, fi parameter space 
(defined in Eq. |TJ for N = 2. Two choices of fi, denoted by 
(a) and (b), were studied in this paper. 



invariance enables us to simulate large lattices and reach 
the close vicinity of the quantum critical point. 

Our key result is that in the broken symmetry phase, 
the scalar susceptibility collapses onto a universal line- 
shape. Its low frequency rise, and its peak, scale together 
with the vanishing Higgs mass mn, which is plotted in 
Fig. [T] The implication of our calculations, is that the 
Higgs mode remains a well defined collective mode all the 
way to the critical point. 

We also determine the amplitude ratios of the Higgs 
mass (in the ordered phase) to the gap A (in the quan- 
tum disordered phase) for the 0(2) and 0(3) models, 
as m H /A = 1.9(1) and m H /A = 2.05(10) respectively. 
These ratios are universal quantities that can be directly 
compared with experiment. Their value differs from the 
mean-field result m^y/A = \ 2, which describes quantum 
critical points in d = 3 for all values of N [19] . 

Model - We consider an O(N) symmetric lattice model, 
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with partition function Z = J Dcj) e Se , where 



Se = 



£n 2 +£(n 2 ) 2 



(i) 



Here, fa is a TV-component real field residing on sites 
i — (x, y, t) of a cubic lattice in discrete Euclidean space- 
time. The model undergoes a quantum phase transition 
at g = g c (p)- See inset of Fig. 1. At weak coupling 
&9 = (.9 — 9c) 1 9c < there is long-range order, and the 
fluctuations include N— 1 gapless Goldstone modes trans- 
verse to the broken symmetry direction {<p). At strong 
coupling Sg > there is a disordered phase with a gap A 
to all excitations. Near the QCP, the long wave length 
properties of this model are captured by an O(N) sym- 
metric relativistic <p field theory in 2+1 dimensions |19j . 
For N = 2 the ordered and disordered phases describe 
the superfluid and Mott insulator of lattice bosons at 
commensurate filling, respectively [2D]. For N = 3, they 
describe the Neel ordered and the gapped singlet phase, 
respectively [T9" l 121 ] I22| . 

Our main focus is the zero-momentum scalar correla- 
tion function in imaginary time, 

Xs (r) = ^E((IW)l 2 l^l 2 )-d^l 2 ) 2 )'( 2 ) 

Xsi^rn) = ^E e ^'" T ^( r )' ( 3 ) 

T 

and the real frequency dynamical susceptibility, given by 



(4) 



Scaling arguments indicate that, near g c , the suscepti- 
bility at small frequencies is of the form |17j 



(5) 



Here, A ~ ,B+|<5g| 1/ is the gap in the disordered phase, v 
is the correlation length critical exponent, and ($+) 
is a universal function of uj/ A in the ordered (disordered) 
side of the transition. The constant C is real, and is a 
regular function of g across the transition. The presence 
of Goldstone modes renders the ordered phase gapless. In 
order to provide a well-defined energy scale that charac- 
terizes fluctuations on the ordered phase Sg < 0, we use 
the gap at the mirror point — Sg across the transition. 
Our goal is to compute the universal scaling functions 
&± and to extract a set of universal parameters that can 
be compared with experiment. 

Methods - The partition function in Eq. ([I]) is refor- 
mulated as a dual loop model, for the cases of N = 2 and 
N = 3 [231. The sum is sampled by a Monte Carlo algo- 
rithm, using the efficient "worm algorithm" [23]. An en- 
hanced performance is achieved by cluster loop updates 
[25j . The simulations were performed on cubic L x L x L 



lattices of size up to L — 200. This allowed us to ap- 
proach within the neighborhood of Sg > 0.36 x 10~ 2 of 
the QCP. 

Our first task was to determine, for each set of pa- 
rameters, the critical point g c . To this end we evaluated 
the helicity modulus p s , as measured by the second mo- 
ment of the winding number in periodic boundary con- 
ditions [24] [26]. At g c , lim£_ i . 00 Lp s (g c ) approaches a 
universal value. g c is then accurately determined from 
the crossing point of Lp s for a sequence of L values [25] . 

The correlation length exponent v is known from pre- 
vious large scale simulations [27] [25]: v-i = 0.6723(3) and 
u 3 = 0.710(2) for N = 2 and N = 3 respectively. 

For N = 2, we study two sets of parameters p\ = 0.5, 
for which g c = 2.568(1), and pi = 2, for which g c — 
3.908(2), see inset of Fig. 1. Thus, the first set of param- 
eters describes a "softer" spin model than the second. 
For N = 3, we use p = 0.5, for which g c = 1.912(1). 

Results - We first extract the gap A in the disordered 
phase at Sg > 0. There, the imaginary part of the scalar 
susceptibility has a threshold at 2 A of the form [T7] , 

$ + (u;/A)~ 9 >/A)e(a;-2A), (6) 

where g 1 {x) — tt/ (in 2 + 7r 2 ^ has a weak logarithmic 
singularity and 7 is a universal constant which equals 1 
in the N = 00 limit. 

The Laplace transform of Eq. |5]) yields the asymp- 
totics of large r as 



X s(t)~A +9 ,(tA)A 4 - 2 ^- 



,-2tA 



tA 



(7) 



where g 1 (tA) is the Laplace transform of g 1 {u>/ A). In 
order to obtain A we fit the simulation results for Xs( T ) 
to Eq. Q. The fit is dominated by the exponentially 
decaying factor, and is insensitive to the parameter 7. 
The critical behavior of the gap near the QCP agrees 
with the form A = B+(5g) u , as shown in Fig. [I] From 
this procedure we extract the values of A + and B + . As 
a check, we also obtained A from the large r decay of 
the single-particle Green's function G(k = 0, r) [2S] and 
found agreement with these results. 

The scaling form of Eq. ^ applies also to Matsub- 
ara frequency correlations x s (iuj m ). Figure [2j a) shows 
the N = 2 scalar susceptibility x s (iuj m ) as a function of 
Matsubara frequency. The data collapses into two curves 
corresponding to $±(zw m /A) on both sides of the tran- 
sition. Collapsing the curves is done by rescaling the 
Matsubara frequencies by A(Sg), fitting the overall con- 
stant shift as a polynomial of Sg, and then rescaling by 
A+A 3 - 2 ^. The black curves in Fig. 2(b) are the rescaled 
functions, which shows good convergence of the numeri- 
cal data in the critical regime. Similar results were ob- 
tained for N — 3 and will be presented elsewhere [23], 

The universality of $± is tested in Fig. [2^b) . The col- 
lapsed functions match closely, especially at small Mat- 
subara frequencies. We use the values of A+ and B+ 
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extracted in an earlier step, without free fitting param- 
eters. This provides a stringent test for the universality 
of the scaling function. 

Analytic continuation - In order to perform the ana- 
lytic continuation of the numerical data, and determine 
the spectral function X^(w) — ImXs(w) we must invert 
the relation 

1 2oj 

Unfortunately the kernel of this integral equation is ill 
posed, which renders the inversion sensitive to inevitable 
numerical noise in x s (iuj m ). To tackle the problem we 
employ the MaxEnt Method [3U] . In this method the in- 
version kernel is regularized by introducing an "entropy" 
functional which is extremized along with the goodness 
of fit. To ensure the validity of the results we track the 
convergence of the MaxEnt spectrum as the statistical 
error decrease for long Monte Carlo simulations. 

Figure [3] presents the spectral function m the 

ordered phase for N = 2 and N = 3. The spectral func- 
tion displays a narrow low energy peak, which softens 
upon approach to the critical point, and a broad high 
energy peak, which does not. This double-peak struc- 
ture is in agreement with the findings of Ref . [TB . The 
position win of the low energy peak as a function of Sg 
is shown in Fig. [l] We find an excellent agreement with 
the expected scaling for the Higgs mode m H = B_\8g\ v 
presented in the red curve. From this fit we extract £>_. 
From the ratio B+/B- we extract the universal ratio of 
the energy scales run /A on both sides of the transition. 

The N — 2 results shown in Fig. grange from \Sg\ = 
0.0036 to \Sg\ — 0.12, corresponding to almost a decade 
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FIG. 2: (a) The scalar susceptibility as a function of Matsub- 
ara frequency for N = 2 and fj, — 0.5. Lower curves, different 
Sg > 0, (disordered phase). Upper curves Sg < (ordered 
phase) (b) Collapse of the same susceptibility onto scaling 
functions $±. Results include two values of {i (see inset of 
Fig.l) 




FIG. 3: Scalar susceptibility x"( w ) (blue online), in arbitrary 
units, for different values of Sg < at /1 = 0.5. Curves shifted 
vertically by Sg. The Higgs energy ran as a function of Sg 
is extracted from the peak positions (dots) as the solid curve 
(red online). 



and a half variation. For the smallest value of 8g, the 
correlation length £ at the mirror point in the disordered 
phase is 22.5 lattice sites. This value satisfy 1 < ^ < L, 
indicating that we are both in the continuum limit and 
in the thermodynamic limit. We have explicitly checked 
that our results do not change upon increasing L. 

Not only does the peak position scale, but the full low 
energy functional form does, as shown in Fig. [4] both for 
N = 2 and N = 3. There we rescale the frequency axis 
by A and the spectral function by A 3 ~ 2 / u to match the 
predicted scaling form of Eq. ^ . Note that the rescaling 
is done without any free fitting parameters since the real 
constant in Eq. ^ drops out from the spectral function. 
The high frequency tail deviates from scaling. This might 
be attributed to contamination from the non universal 
part of spectrum and systematic errors introduced from 
the MaxEnt analysis. The observed functional scaling 
demonstrates that the Higgs peak is a universal feature 
in the spectral function that survives as a well-defined 
excitation arbitrarily close to the critical point. 

The peak position in units of A is shifted to higher en- 
ergies for the N = 3 case compared to N = 2. This trend 
agrees with the prediction made in |17j that m#/A in- 
creases monotonically with N. We also obtain the fidelity 
F = m^y/r, where T is the full width at half-maximum. 
We measure T with respect to the leading edge at low 
frequency, since at low frequencies there is less contam- 
ination from the high frequency non-universal spectral 
weight. Since the entire functional form of the line shape 
is universal, F is a universal constant that characterizes 
the shape of the peak. We find F = 2.6(1) for N = 2 
and for F = 2.25(10) for N = 3. 

In the ordered phase, the asymptotic low frequency rise 
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of the susceptibility was predicted [TH [T51 E] to be 

~ (lu/A) 3 , u«A«l, (9) 

The oj 3 rise is due to the decay of a Higgs mode into a 
pair of Goldstone modes. On the other hand, Fig.[4]does 
not display a clear a; 3 low frequency tail. An alterna- 
tive method to look for this tail exists, without the need 
to analytically continue the numerical data to real time. 
Equation ^ transforms into the large imaginary time 
asymptotics Xs (t) ~ 1/t 4 - 

For N = 3 we confirm the asymptotic behavior of 1/r 4 
within the numerical limitations. Interestingly, for N = 2 
we do not find a conclusive asymptotic fall-off as 1/t 4 , 
Instead, the data fits better to an exponential decay, as 
in the disordered phase (see Eq. Q), although the power 
law might have a small amplitude below our statistical 
errors. In both cases, we can safely conclude that the 
spectral weight of the Higgs peak dominates over the low 
frequency w 3 tail, enhancing its visibility. The large r 
analysis is discussed elsewhere [25] . 

Discussion and Summary- Our results are directly ap- 
plicable to all experimental probes that couple to a func- 
tion of the order parameter magnitude. For example, the 
lattice potential amplitude in the trapped bosons sys- 
tem [IB] [IS], or pump-probe spectroscopy in Charge Den- 
sity Wave systems 6-8 . Such a probe can be expanded 
near criticality in terms of the order parameter fields and 
their derivatives, 

G(x, r) = a\$\ 2 + (3\dJ\ 2 + 7 (|^| 2 ) 2 + . . . (10) 

So long as a ^ 0, the first term is more relevant than the 
rest. Hence, the scalar susceptibility defined in Eq. ^ 
dominates the experimental response at low frequencies 
and wave vectors. 

In summary, we have calculated the scalar susceptibil- 
ity for relativistic 0{2) and 0(3) models in 2+1 dimen- 
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FIG. 4: Rescaled spectral function vs. cj/A for N = 2, 3, 
at u — 0.5. At low values of oj/A, these curves collapse to 
the universal scaling function $" (u/A), in accordance with 
Eq.[5] 



sions near criticality. We have demonstrated that the 
Higgs mode appears as a universal spectral feature sur- 
viving all the way to the quantum critical point. Since 
this is a strongly coupling fixed point, the existence of 
a well defined mode that is not protected by symmetry 
is an interesting, not obvious, result. We presented new 
universal quantities to be compared with experimental 
results. 
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